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Abstract 

We describe a limiting case when nonlinear dynamics of 
an ac-driven semiconductor superlattice in the miniband 
transport regime is governed by a periodically forced and 
damped pendulum. We find analytically the conditions 
for a transition to chaos and consider an influence of 
temperature on the effect. We also discuss fractional dc 
voltage states in a superlattice originating from phase- 
locked states of the pendulum. 

PACS: 73.21.Cd; 72.20.Ht; 05.45.-a; 05.45.Ac; 
72.30.+q 

Keywords: Semiconductor superlattice; 
chaos; phase-locking; Josephson junction 



1 Introduction 

Strongly nonlinear dynamical effects in miniband semi- 
conductor superlattices (SSLs) driven by an intense high- 
frequency electric field attract growing interest in the last 
few years @, §, § |, g, |, @, |, §. This theoretical ac- 
tivity, which includes investigation of dissipative chaos 
pi [/J and a spontaneous generation of dc voltage in 
a purely ac-driven SSL |, | |, |, @, |, |, |o| , has been 
stimulated by a recent progress in the experimental stud- 
ies of nonlinear electron transport in the SSLs driven by 
THz field [[n]. On other hand, recent experiments fl2|| 
reveal an important role of coherent plasma oscillations 
in nonlinear dynamics of carriers in SSLs. 

The counterparts of the effects of chaos and sponta- 
neously dc voltage generation in the THz-driven SSL are 
known in the rf-driven Josephson junctions for a rather 
long time |l| ||, [l| ||, (for a review, see @). So- 
called an inverse ac Josephson effect (a quantized dc volt- 
age generation) [TtJ has already found applications 
in the design of a zero-biased Josephson voltage stan- 
dard [ fi"8| , which offers several advantages in comparison 
with the conventional voltage standards based on an ac 
Josephson effect 
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A transition to dissipative chaos 
often limits the stable operation of such kind of voltage 
standard . Nonlinear dynamics of driven plasma os- 
cillations in a junction is well-described by the fairly pop- 
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ular Resistively Shunted Junction (RSJ) model |2(| |21| ], 
which is equivalent to a periodically forced and damped 
pendulum - one of the most important paradigms in the 
modern nonlinear science p2]| . 

Within the semiclassical approach [^3|, ac-driven 
plasma oscillations in a single miniband SSL are de- 
scribed by the set of balance equations for the average 
electron velocity (current), the average electron energy 
and the electric field (voltage) generated by an electron 
motion in SSL g, g, It has been shown in 

rcf. H that this set of three nonautonomous nonlinear 
first-order differential equations is isomorphic to the pe- 
riodically forced Lorenz model 24 . In comparison with 
the RSJ model, this system is more complex from a view- 
point of nonlinear dynamics and until recently it has 
been analysed mainly numerically |l] [| |[ f|, 0, [|. 

However, it should be noticed that a pendulum-like be- 
havior of miniband electrons in SSL has been established 
earlier for some limiting cases. First, plasma oscillations 
in SSL have been described within a collisionless approx- 
imation by a pendulum without damping in p5| , |26| , |^| . 
Second, describing the frequency-to-voltage conversion 
in ac-driven SSL, Dunlap et al. 27| have found that 
their model of SSL demonstrates a behavior similar to 
the phase-locked solutions of damped pendulum Jl(| for 
u > ujpi (oj and u) p \ are the ac electric field frequency 
and the miniband plasma frequency, correspondently). 
These earlier findings indicate that dynamics governed 
by a pendulum can be underlying for nonlinear dynamics 
of ac-driven SSL. 

In this work, we show that the driven damped pendu- 
lum could describe a stationary nonlinear dynamics of 
an ac-driven SSL in the limit when a characteristic scat- 
tering constant of electrons with impurities and phonons, 
7, is less than an an ac frequency u>. In contrast to the 
RSJ model, where the voltage across junction is propor- 
tional to the velocity of pendulum, the voltage across 
SSL depends on both the coordinate and the velocity 
of pendulum. These findings possess us to apply several 
exact and approximate results obtained earlier for a pen- 
dulum to an analytic investigation of regular and chaotic 
dynamics in SSL. Namely, we will show the following: 

(1) Chaotic dynamics is natural for a relatively low 
damping, i.e. for 7/wpi <C 1. 

(2) The SSLs with wide minibands are more stable 
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against perturbations leading to a transition to chaos 
in comparison to the case of narrow miniband SSLs. 

(3) The Mclnikov criterion of transition to chaos in the 
pendulum qualitatively explains the position of low fre- 
quency boundary of transition to chaos in SSL. 

(4) For a relatively low damping, 7/cjpi <C 1, phase- 
locking in driven pendulum results in a generation of dc 
bias across of SSL; the dc voltage per period of SSL is 
approximately proportional to a half-integer multiplier 
of fku. 

(5) The difference between generated dc voltage per su- 
perlattice's period and hw is proportional to (7/cjpi) 2 at 
7/wp] < 1. 

(6) For a fairly strong scattering, phase-locking in driven 
pendulum may not result in a generation of quantized (or 
an approximately quantized) dc bias across SSL. 

Our present analytic research explains and further de- 
velops the results obtained numerically in the previous 
publications || | |, §. 

2 Pendulum limit in superlattice 
balance equations 

We consider electron transport through a single mini- 
band of a spatially homogeneous SSL under the influ- 
ence of ac current. We assume the tight-binding depen- 
dence e{k) = A/2[l — cos(fca)] of the electron energy e 
on its quasimomentum k along SSL's axis, where A is 
the miniband width, and a is the period of SSL. The 
electron dynamics could be described by the following 
balance equations |], ||, @, H [J (we use notations of 
ref. §) 



The electron energy within miniband relaxes to the ther- 



v 
u 



uw — 7W, 

—uv — j(w — w 

lUqV + Iq sinwi, 



(1) 



Here the scaled variable u = eaE(t)/h describes the elec- 
tric field E (or the voltage per SSL's period), generated 
by an electron motion along the SSL axis; v — m^Va/h 
is the scaled electron velocity (V is the electron velocity 
averaged over the time-dependent distribution function 
satisfying the Boltzmann equation, mo = is the ef- 
fective electron mass at the bottom of miniband) and 
w = (e — A/2)(A/2) _1 , where e is the average elec- 
tron energy within the miniband. Following this scaling 
for the energy, the lower (upper) edge of the miniband 
corresponds to w — —1 (w = +1). Next, u) is the fre- 
quency and Jo is the scaled amplitude of driving ac cur- 
rent; 7 is the phenomenological relaxation constant de- 
scribing different channels of electron scattering in SSL; 

1 /2 

ujq = uj p \ (C s i/C) ' is the generalized plasma frequency, 

1/2 

where uj p i — (47re 2 iV/moeo) is the miniband plasma 
frequency, N is the number of electrons per unit volume, 
eo is an average dielectric constant for the SSL material, 
C s i = ^4 is the self-capacitance of SSL with the length 
I and the cross-section area S, and C is a parasitic ca- 
pacitance, which is in parallel to the SSL's capacitance. 



mal equilibrium energy value, w. 
of a nondegenerate electron gas 
temperature dependence |M Blf 



'eq (w eq < 0); in the case 
u cq has the following 
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where Io,i are the modified Bessel functions, T is the 
lattice temperature, and ks is the Boltzmann's con- 
stant. Note that the dependence (||) well describes 
temperature-induced modifications of the Drude conduc- 
tivity in the narrow miniband SSLs as have been ob- 
served in the experiments [^8[ ^9) . 

First two equations of the set ([l]) are the well-known 
balance equations of Ignatov and Romanov 
which could be derived from the Boltzmann transport 
equation with a constant relaxation time 7. Thus, 
instead of finding the time-dependence of distribution 
function f(p,t) (p = Hk) from the Boltzmann equation 
and further calculation of the average electron velocity 
V = J dpf(p, t)de{p) I dp and the average electron energy 
e = J dpf(p,t)e(p), it is possible just to solve the cor- 
responding balance equations for the average variables. 
Both these approaches give the same result for the tight- 
binding dispersion relation. 

The third equation in (Q) could be obtained from the 
Maxwell or Kirchoff equations §, || §, @, §]. To the 
best our knowledge, three coupled balance equations in 
the form ([!]) have appeared first in the work of Tetervov 
J32| , who used these equations for a description of decay 
of undriven (Jo = 0) nonlinear plasma oscillations in 
SSL. 

Now we want to show how the pendulum limit arises 
in the superlattice balance equations (|l|). We introduce 
new variables A and 9 as w — —AcosO, v = —As'm9. 
We can re-write Eqs. ([l]) in the form (see Appendix |A|) 

6 + Au$ sin 9 + (6, A)6 = 

I Q smut + ~/ 2 w cq A- 1 F(9,A), (3) 
A = — -fA — 7w cq cos 6*, (4) 



where 



G{9,A) = -WeqA' 1 cos6, (5) 
F(9,A) = sm9 + w eq A~ 1 sm6cos9, (6) 
9(t) = u(i) + 7™ oq A- 1 (i)sin(9(i). (7) 

Formally, the set of equations (0)-(0) is not simpler than 
the original balance equations (Q). However, this form is 
more suitable for the consideration of pendulum limits. 

We start our analysis with the case of collisionless SSL, 
7 = 0. As follows from Eq. (f|), A 2 = v 2 (t) + w 2 (t) is 
the integral of motion in this case. For natural initial 
conditions v(0) — 0, w(0) — w cq , we have A = —w cq . 
The balance equations are reduced to single equation of 
undamped and periodically forced pendulum 



* + (- 



(; eq Vo sm ^ = -^0 sin u)t, 



(8) 
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where v(t) — w cq sin9(t), w(t) — w cq cos9(t), u(t) = 
8{t). Earlier the same equation has been obtained within 
the semiclassical Hamiltonian approach by several au- 
thors (25[ |26|]. It was used in the investigations of non- 
linear plasma oscillations |2jJ, the instabilities and the 
Hamiltonian chaos p6| in ideal ac-driven SSLs. 

Eq. (||) is valid only for a short time t <C 7 _1 - Now we 
want to show that a pendulum equation can also describe 
electron's dynamics in the ac-driven SSL at long time 
(f 3> 7 _1 ) providing lo ^> 7. Eq. represents a pen- 
dulum, parameters of which are varied self-consistently 
due to their dependences on A and 9. To begin with 
we limit our attention to the time interval T = 2tt/lo. 
During this characteristic period of oscillations (or rota- 
tions), T, the variable A only lightly decays for lo 3> 7 
(see Eq. (|)). Now we compare characteristic damping 
terms in Eqs. d) and (|). Following Eqs. (g) and (§), 
the averaged over period of ac field dissipative term is 



(^G(6,A)) = - 1 w cq {A- 



cos ( 



Thus, the ratio of characteristic dissipative terms in 
Eqs. (g) and in Eq. (Q), calculated during the time in- 
terval T, is following 



( 7 G(9,A)) 
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w eq (A~ 



cost 



< 1 



This ratio is small because the term cos 9 changes its 
signs (oscillates), while A stays almost constant. Next, 
we consider more long time interval of more number of 
periods T and come to a conclusion that a dissipation 
in Eq. (^) is, in effect, much more strong than a dissipa- 
tion in Eq. (||); the variable A(t) approach its stationary 
value faster than the variable 9(t) approaches its attrac- 
tor. Therefore, if additionally the amplitude of driving 
current is not large, Iq/cOq < 1, the variable A can be adi- 
abatically eliminated p3) . That is, we can substitute the 
stationary value A(t) — > A Et (t) — —w cq cos9 in Eq. (|3|). 
Observing that G(9, A st ) = 1 and F{0, A st ) = 0, we find 
that Eqs. (g), (@) take the form 

tp + jtp + (—w cq )u)Q sin-0 = 2Jq sin tot, (9) 



u = vV2 + 7tan(V>/2), (10) 

where we introduced ip(t) = 29(t). This pendulum rep- 
resentation of superlattice balance equations constitutes 
the main result of this paper. 

It is useful to compare some well-known solution of 
Boltzmann equation for the tight-binding lattice and the 
corresponding solution of pendulum equation (g). This 
test can be simplified employing an analogy with the RS J 
model of Josephson junctions, which is also governed by 
a periodically driven and damped pendulum. Consider 
a SSL driven by a given ac electric field E — Eq cos cut 
with the frequency being much greater than the mini- 
band plasma frequency, w > wj. In this case we have 
Iq = lo s lo with lo s — eE a/h [g, f7| || ||. For such I Q and 
lo, we find that the solution of pendulum equation for the 



variable v(t) = sin 9 cos 9 is same as the well-known ex- 
pression for the time-dependent electron's velocity v(t) 
derived from the exact solution of the Boltzmann equa- 
tion |3(], if lo 3> 7 and io/^o < 1 ( see Appendix 
[b]). Therefore, we see that the pendulum representation, 
Eqs. (p|)-(|l0|), gives same result as the Boltzmann equa- 
tion within the range of its validity. Moreover, this range 
of parameters, Iq/uq < 1 and to loq, corresponds to 
z = lu s /lo <C 1 for w>wo, what are quite typical values 
of the field strengths and the frequencies in the modern 
experiments on harmonics generation and detection of 
THz radiation in SSLs |ll[] . Finally, it is worth to notice 
that from the viewpoint of analogy with the Josephson 
junctions, a described generation of harmonics of an ac 
electric field in SSL is completely equivalent to the stan- 
dard ac Josephson effect jl8[ pif ) i.e. to a generation 
of current's harmonics under the action of a sinosuidal 
voltage. 

In this paper we are mainly focused on nonlinear dy- 
namical effects arising within the pendulum limit. With 
the reference to the pendulum representation of balance 
equations, we can make the important conclusion: The 
thresholds of instabilities in an ac-driven semiconductor 
superlattices, including the transition to chaos, are same 
as in the RSJ model of an ac-driven Josephson junction 
providing lo 3> 7. 

Before we proceed with a consideration of instabilities 
and chaos, we should note that the SSL balance equa- 
tions (||)-(|7]) can be also simplified in the opposite limit- 
ing case of frequent collisions, to <C 7. Really, assuming 
\8\ ~ to < 7 we see from Eq. (||) that A(t) follows 9(t) 
adiabatically. Therefore Eqs. (§), (§), (§), (0) take the 
form 

'ftp + loq sin tp = 2I sin cot (11) 

with u = 7tan(?/>/2). Such kind of the first order equa- 
tion is well-known in the theory of Josephson junctions 
p5[ pl| . Importantly, it also arises in the description of 
an interaction of electromagnetic radiation with the lat- 
eral SSLs (l0[ [36|. Chaos is impossible in the first order 
differential equation p2] |. However, the overdamped pen- 
dulum ( pd| ) driven by a two- frequency quasiperiodic force 
still can demonstrate very complex dynamics known as 
a strange nonchaotic attractor j37|. 

Now we turn to the consideration of conditions for a 
transition to chaos and an influence of temperature ef- 
fects using the pendulum representation of balance equa- 
tions @, ©. 

3 Periodicity and chaos 

The theorem of Levi (38) states that the stationary dy- 
namics of a periodically forced pendulum is never chaotic 
in the overdamped limit, which in our case reads 



C = 



7 



11/2 



> 2. 



L0 



(12) 



The numerical simulations of the driven damped pendu- 
lum performed for the wide ranges of Iq and lo reveal the 
absence of chaos already for C ~ 1 pgt] . 
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For a narrow miniband at room temperature, the ther- 
mal equilibrium energy becomes close to the center of 
miniband, i.e. \w cq \ is small (see Eq. (^)). Therefore, 
the criterion for absence of chaos can be easier sat- 
isfied in the SSLs with narrow miniband in a compari- 
son to the case of wide miniband SSLs. In particular, 
for ksT » A we have w cq « — A/2fcsT and there- 
fore C oc 7/ A (really, coo oc m^ 2 oc A 1 / 2 providing 
ItUeql 1 / 2 ^ oc A). Thus, with a decrease of miniband 
width A, the value of C oc 7/A increases making the 
condition ( |l2| ) to be more robust. 

As an example we consider the narrow miniband SSL 

with A = 3 meV, a = 22 nm, N = 2 x 10 15 cur 3 [§§] 

and 7 _1 w 0.5 ps E9], the miniband plasma frequency 

is w p i = 2.1 x 10 1 rad/sec and 7/u/pi « 1. For the 

temperature T = 50 K, w cq « —0.33 and the factor 
1/2 

w q I is 1.7. For the higher temperature, T = 200 K, 
-1/2, 



3.5 and chaos is always impossible. However, 
chaos is still possible for a wider miniband and a longer 
relaxation time (see corresponding estimates in Refs. 0, 

The criterion of absence of strange attractor, Eq. (112), 
has been obtained within the pendulum limit and thus 
is formally valid only for w>7. However, summing up 
the results of numerical simulations of the superlattice 
balance equations, we speculate that strange attractor 
does not exist at any frequency w, if the inequality (|l2| ) 
is satisfied. 

To find an analytic criterion of a transition to chaos in 
the pendulum limit, we apply the method of Melnikov 
[ [io[ . This method has been used by several authors for a 
determination of the conditions of transition to chaos in 
the periodically driven pendulum (or for the RS J model) 
[ ^5| fL8| . In our case, the Melnikov's criterion for a tran- 
sition to chaos is 



Io > I c Q r 



27^0 Kq 1 1/2 



cosh 



2 \w eq \V 2 u 



(13) 



(cf. Eq. (60) in Ref. JlJ]). This formula is derived with 
the assumption that conditions 
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«1, T 1I^ 1 



|lO e q| 1/2 Wo koqlw 2 

are satisfied. For u> < luq Eq. ( p"3| ) takes the form 
2 



'-juj \w cq \ 1/2 , 



(14) 



In the case of SSL driven by the ac electric field E — 
Eq coswt, one can get Iq — lj s uj with uj s — eE^ajh |2|, |?], 
||, S|. Then, the formula ( fil| ) gives for the boundary of 
chaotic region in the uj s -lu plane in the form 



1 1 1 1/2 

UJ S UJ = -7W |Weq| ' . 



(15) 



Earlier, the low-frequency boundary of chaos in the sim- 
ple form uj s uj — const has been found numerically in the 
works j^, ||, [ll| , but has been unexplained until now. Eq. 
L5|) can qualitatively explain this dependence. 



4 Phase-locking in pendulum and 
dc voltage generation in super- 
lattice 

Phase-locking is another nonlinear dynamic phenomenon 
which an ac-driven pendulum can demonstrate along 
with chaos [Q. In the physics of Josephson junctions, 
phase-locking at zero bias is known as the inverse ac 
Josephson effect Jll], |l7], [l8| Here we are dis- 

cussing the influence of phase-locking in the pendulum 
on the physical properties of SSL. In the driven pendu- 
lum phase-locking means that 



(16) 



where n and I are integer numbers and (...) stands for 
the time averaging over the period of ac drive 2tt/uj 
|i~4|, ^8). Majority of stable phase- locked states are inte- 
ger (n ^ 0, I = 1); however, a fractional phase-locking 
(I > 1) can also exist. In the case of Josephson junc- 
tion, the voltage U across the junction is proportional 
to ip; therefore phase- locking results in a generation of 
quantized dc voltage (U) oc (n/l)u 0, |l|, Q . 

For the superlattice problem, the voltage is a function 
of both the velocity and the phase of the pendulum (see 
Eq. ([lOl)). Therefore, phase-locking in the pendulum 
( |16|) determines a generation of dc voltage in the SSL as 



Rnl 



nut 7 „ 

1+ J-J^ 
/ 2 \ o^o 



(17) 
(18) 



For a low damping 7/wpi 1, we immedeatly get from 
( p"7j ) the dc voltage {u) « (n/2l)uj. Such kind of almost 
half-integer dc voltage states in a pure ac-driven SSL 
have been recently observed in the numerical simulations 
of balance equations [Q [)| . 

The formula ( |l7j) tells us nothing about a stability of 
dc voltage states in SSL that correspond to the phase- 
locked states of the pendulum with different n and /. 
However, numerical simulations demonstrate that the dc 
voltage states, which are close to the integer states (i.e. 
to the states with I = 1 and n being even integer), are 
more typical then fractional states || [1. 

What is important, the dc voltage per superlattice pe- 
riod is not exactly (n/2l)to even for a weak damping. 
Moreover, a weak dependences of (u) on the ac ampli- 
tude |Q, |d| and the frequency |l] |4l]] have been found 
numerically at 7/wpi <C 1. The appearance of such kind 
of dependences within the pendulum limit can be under- 
stood, if one takes into an account that (u) is a function 
of both the time-average velocity and of the time-average 
coordinate of pendulum (see Eq. (p"7|)). Phase-locked ro- 
tational state of the pendulum can be represented in the 
form [[J 



■0(f) = ipo + noj + a p sm(pujt + fi p ), (19) 
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where V'Oj ol p and /i p are constants and we restrict our 
consideration to the leading phase-locked states, (ip) = 
nuj. The amplitudes a p and the phases fi p in Eq. (|l^) are 
functions of the main physical parameters of our prob- 
lem, i.e. Iq, ui and 7. Substituting the expression ( |l9| ) in 
Eq. (18) one can see that the dependence of dc voltage 
(u) on the ac amplitude and frequency arise in Eq. 
via the term R n i. 

It is important to know how this "dissipative cor- 
rection" to the quantized dc voltage, {u) — nuj/2 = 
(7 /too)R n , scales with an increase of dissipation. To es- 
timate the dependence of R n on 7/wo for several lowest 
phase-locked states of the pendulum, we consider the 
limit of fast rotations a< 1. In this case formula ( |l9| ) 
takes the form 



?/>(£) = tpo + nuj + a smicot + /i). 



(20) 



Substituting Eq. ( |20| ) in Eq. (^|) and equating zero har- 
monics, we have the formula (cf. (Ml (42[) 



Wq J n (ct) sin(n^) = —n^/uj, 



(21) 



which gives us an opportunity to understand the char- 
acteristic dependences of a and ip on 7 and uj. Next, 
substituting Eq. (|2(]) in Eq. ( |l8| ) and calculating inte- 
grals in the limit a C 1, we obtain (see Appendix ^]) 

(tan (ip/2)) n=1 = asm(ri + 0(a 2 ), (22) 
(tan(^/2)) n=2 = -0.25a 2 sin(2^) + 0{a i ). 



Now combining Eqs. (|18j), (|21|), ( |22[ ) and using the 
asymptotics of the Bessel functions J n (z) ks z"/(2"n!) 
(z <C 1), we find 



Rn - (-1) 



n f7 
W ' 



(n=l,2) 



(23) 



Thus, for a weak dissipation (7/wpi <C 1), the dc voltage 
spontaneously generated in a pure ac-driven SSL is 



<u) = iuj/2 + 0( 7 2 /^pi) 



(24) 



and it has a quadratic dependence on 7. 

With a further increase of 7, the dissipative term in 
Eq. ([IT]) becomes more important. Moreover, for a large 
dissipation the phase-locked rotational states (|l9|), ( |20| ) 
do not exist in the pendulum anymore fl^] . However, the 
unquantized dc voltage still can be generated in SSL due 
to different mechanisms, as has been demonstrated nu- 
merically in H . The discussion of relationship between 
a pendulum dynamics and a dc voltage generation in a 
superlattice for an arbitrary damping will be presented 
elsewhere. 



5 Discussion and conclusion 

It is instructive to compare the pendulum equation 
derived from the superlattice balance equations in the 
limit of rare collisions (7 <C uj) and corresponding motion 
equation that can be obtained employing the Newton 



law. We introduce simplest form of phenomenological 
dissipative term to the motion equation for electron's 
momentum as 

P + IP = —eE(t). 

Further, combining this equation and the standard def- 
inition of electron's velocity V — e(p)/dp, the tight- 
binding dispersion relation, and the Maxwell equation 
(third equation in the set (Q)), we get the pendulum 
equation in the form 

(f> + 70 + ujq sm = ^0 sin cut, 



where <f> = pa/h, \ 



II' c 



1 for simplicity. In this ap- 



proach the scaled voltage is u = <f). Certainly, this equa- 
tion is similar to Eq. (^|). The difference arises mainly 
in the dependence of voltage on coordinate and momen- 
tum of the pendulum (cf. Eq. (^0|)). In particular, 
the voltage is dependent not only on the momentum of 
pendulum but also on its coordinate in the Boltzmann 
equation approach. 

The crude approach based on the Newton law uses 
a trajectory of an individual electron. Therefore, this 
approach can not take into an account the changes in 
a distribution function of electrons induced by an elec- 
tric field. In contrast, the approach based on the Boltz- 
mann equation does include such kind of corrections in 
the term proportional to 7. In the limit of infrequent 
collisions, 7^0;, this term can be often ignored be- 
cause 9 ~ to S> 7. However, the dissipative contribution 
to voltage can become important if one considers a dc 
voltage generation, (u) ^ 0. As we have demonstrated in 
section ^, this dissipative correction can be important in 
the consideration of phase locking, (0) = uj. The contri- 
bution of dissipative term to the dc voltage, (w), can be 
even more important in the case of lateral superlattices, 
for which (6) =0 10. 

In summary, we show that pendulum equation can 
well describe the dynamics of electrons in a miniband 
of ac-driven semiconductor superlattice in the limit of 
weak scattering. We demonstrate that several earlier 
numerical findings concerning chaos and spontaneous dc 
voltage generation in the superlattices can be explained 
within such approach. We establish a link between the 
theory of diffusive transport in ac-driven superlattices 
and the theory rf-driven Josephson junctions. We point 
out that dissipative terms arising in the dependence of 
voltage across superlattice on the momentum and the 
coordinate of the pendulum should play an important 
role in many situations. 

We believe that described in this work simple and in- 
structive theory can be useful in an explanation of the 
results of experiments |ll], |l2) as well as that it can stim- 
ulate other experiments on the THz-field induced non- 
linear phenomena in superlattices. 
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A Appendix. 

In this appendix we will derive the equations (|jj|), (Q), 
© > @i @ • begin with we introduce a complex vari- 
able Z = w + iv and represent the superlattice balance 
equations M) in the form 



Z 



-iuZ — jZ + 7W 



(25) 



Substituting the parametrization Z = —Aexp(i9) (i.e. 
to = — A cos 8 and v = —A sin 8) in (|25| ) and equating 
real and imaginary parts, we get 



A = —jA — 7w eq cos 8, 
8 = u + 7w oq j4 _1 sin 6?. 

From the last equation we have 

8 = u + 7 w cq^ sm #] 
Taking into account that 



(26) 
(27) 



(28) 



dt 



(L4 -1 sin0)] = -A~ 2 Asm6 + A^OcosO (29) 
-fA^ 1 sin 8 + jw cq A~ 2 sin 8 cos 8 + 8A~ X cos 8 



and using the third equation of set ([!]), we immediately 
get from Eq. @ the Eqs. (|), (|), (§). Finally, it is 
easy to see that Eqs. (|2^) and (|27j) are same as Eqs. (|J) 
and (m), correspondently. 



B Appendix. 

In this appendix we will show that well-known result 
of Tsu-Esaki-Ignatov-Romanov from the theory of ac- 
driven SSLs 



34 



can be reproduced within the pen- 
dulum representation of the SSL balance equations (j^), 
©• 

Solving the Boltzmann equation for the tight-binding 
lattice driven by the electric field E = Eg cos u>t with 
to 3> wo > Ignatov and Romanov ]3C| ] have found for )<w 
the following time-dependence of the averaged electron's 
velocity 



term (oc 8 ~ ui). Further, we can employ an analogy with 
the RSJ model of Josephson junctions. It is known for 
the RSJ model that the sinusoidal voltage is the solution 
of pendulum equation if u> 3> ujq and 7 <C u> JT8| . In 
our case it means that 8(t) = zsiriujt is the solution 
of Eq. (j^f. Further, substituting this time dependence 
of 8(t) in the expression for averaged electron's velocity 
v = sin 8 cos 8, we can determine the current generated 
in the SSL under the action of the ac field. (Of course, 
it is the analog of the ac Josephson current ) . Using 
the Bessel expansions 



cos(zsini) = J (z) + 2^_] J2k{z) cos2/ci, 
fc=i 



sin(z sin t — 2 Jik-\ ( z ) cos(2fc — l)t, 
fc=i 

we find the averaged electron's velocity in the form 

00 

v(t) = 2J {z)^2j 2 k-i{z)sm(2k-l)ujt (31) 



fc=i 



+ 2 ^2 J2k 1 ( z ) J 2k 2 -l(z) X 
fei,fc 2 >l 

[sin (2(fc 2 - fci) -l)ut + sin (2(fc 2 + fa) - 1) wi] 



Supposing that 



we have 



C0 S < 



(32) 



Z EE — < « 1 

because w/cjo 3> 1. Now it is easy to see that for z <C 1 
Eq. © is same as Eq. @). 

Thus, we have showed that for 7 <C a; and, addi- 
tionally, if the condition (^) is satisfied, the pendulum 
equation can reproduce the result of Tsu-Esaki-Ignatov- 
Romanov. Finally, combining the definition Iq = lu s lo 
and Eq. (|32|), we see that the inequality ( |3^ ) just gives 
another condition of the validity of the pendulum ap- 
proach: Io/ujq < 1. 



C Appendix. 

In this appendix we will derive the equations 
need to calculate the integrals 



(©■ We 



00 1 f 27r rnx a \ 

;(t)=^a 2fe _ lS in(2fc-l)o;t, a 2k . x = 2J (z)J 2k -i(z), < tan (W2)} = ^ J q dxtan + - sin(x + M )J , 

k=l ( Q' 



(30) 

where z = lo s /lo = (eEoa)/fko. Same result has been 
obtained earlier by Tsu and Esaki for the limiting case 
z«l|. 



(33) 

where n = 1,2. In the limit a < 1 we can use the 
expansions 



Now turn to the pendulum equation. Hence we have 
the ac field of a very high frequency, we can neglect the 
last term (oc 7) in Eq. ( |l0| ) in comparison with the first 



tan 



(I + fria(* + AO)=tan(f) 



sin(a: + fj,) 
cos(x) + 1 



a-O (a 2 ) , 



(34) 



G 



/ a . . A , a sin(x + /i) 

tan x + — sm(x + fi) ) — tan(x) + — 

V 2 / 2 cos ^(x) 

a 2 (— 1 + cos 2 (a; + /i)) sin(x) 

4 cos 3 x 

a 3 (—1 + cos 2 (x + //)) sin(x + //) (2 cos 2 x — 3) 
24 cos 4 x 

+0 (a 4 ) . 



After integration of the expression (t 
tain Eqs. (||). 
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